ON THE DRAWDOWN OF COMPLETELY ASYMMETRIC LEVY PROCESSES 
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Abstract. The drawdown process y of a completely asymmetric Levy process X is equal to X 
reflected at its running supremum X: Y — X — X . In this paper we explicitly express in terms 
of the scale function and the Levy measure of X the law of the sextuple of the first-passage time 
of Y over the level a > 0, the time Gt„ of the last supremum of X prior to Ta, the infimum X ^ 
^ , and supremum Xr^ of X at Ta and the undershoot a — Yt^- and overshoot Yt^ — a of Y at Ta. As 

I application we obtain explicit expressions for the laws of a number of functionals of drawdowns and 

rallies in a completely asymmetric exponential Levy model. 
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1. Introduction 



A completely asymmetric Levy process is a real-valued stochastic process with cadlag paths 
that has independent stationary increments whose jump sizes all have the same sign. Its drawdown 
process, also known as the reflected process, is the difference of its running supremum and its current 
^SJ ' value. Closely related is the rally process which is defined as the difference of the current value and 
the running infimum, and is equal to the drawdown process of the negative of the process. This 



paper is concerned with a distributional study of the drawdown process of a completely asymmetric 
Levy process For such a Levy process X the law is identified of the following sextuple concerning 
. its drawdown process Y: 

o _ _ 

(1.1) {Ta,Gr^,Xr^,X^^,a-Yr,_,Yr,-a), where Yt = sup X^ - Xt, t>0. 

0<s<t 

Here the components of the vector are given by r^, the first-passage time of Y over a level a > 0, 
^ ■ Gr^, the last time that X is at its supremum prior to Tq, Xt^ and 2£.r^, the running supremum and 
infimum at Ta, and a — l^a- and Yr^ — o, the undershoot and the overshoot of Y at the epoch Ta- 

The drawdown process has been the object of considerable interest in various areas of applied 
probability. It has been studied for instance in queueing theory {e.g., Asmussen [5]), risk theory and 
mathematical genetics. The drawdown process has also been employed in financial modelling, in the 
construction of tractable, path-dependent risk/performance measures. In the context of real estate 
portfolio optimisation, Hamelink & Hoesli [20] considered the running maximum of the drawdown 
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process as an investment performance criterion. Checkhlov et al. [T^ introduced a one-parameter 
family of portfolio risk measures that was called conditional drawdown and defined to be equal 
to the mean of a percentage of the worst portfolio drawdowns. Pospisil et al. [32] proposed the 
probability of a drawdown of a given size occurring before a rally of a given size as risk measure, 
and calculated this probability in the setting of one-dimensional diffusion models; the finite horizon 
case for Brownian motion was treated by Zhang &: Hadjiliadis |37| . 

A number of papers has been devoted to a distributional study of functionals of the drawdown 
process. The joint Laplace of the time to a given drawdown and the running maximum of a Brownian 
motion with drift was derived by Taylor [36]; this joint law was obtained by Lehoczky [21] in the case 
of a general diffusion. An explicit expression for the expectation and the density of the maximum 
drawdown of Brownian motion was derived by Douady et al. [16] ; the case of Brownian motion with 
drift was covered by Magdon et al. [26] where also the large time asymptotics of the expectation 
were derived. 

The drawdown process also features in the solution of a number of optimal investment problems. 
Under the geometric Brownian motion model the optimal time to exercise the Russian option, which 
pays out the largest historical value of the stock at the moment of exercise, was shown by Shepp & 
Shiryaev [M] to be given by the first-passage of a drawdown process over a certain constant level. 
Such a first-passage time is also optimal when linear cost is included (Meilijson [27]), or under a 
spectrally negative Levy model for the stock price (Avram et al. [3]). 

Although still widely used as benchmark, mainly on account of its analytical tractability, it is 
by now well established that many features of Samuelson's classical geometric Brownian motion 
model for the price of a stock are not supported by empirical data. A class of tractable models that 
captures typical features of stock returns data such as fat tails, asymmetry and excess kurtosis is 
that of exponential Levy processes. This class has received considerable attention in the literature 
— we refer to Cont & Tankov [15] and Boyarchenko & Levendorskii [8] for background and references. 
By restricting ourselves to Levy processes with jumps of a single sign, we are able to draw on the 
fluctuation theory for this class of stochastic processes, which is considerably more explicit than 
in the case of general Levy processes. Empirical support for a model from this class was given in 
Carr & Wu [9], where options on the S&P 500 index were studied. Carr & Wu [9] demonstrated 
that the finite-moment log-stable model, which is an exponential Levy model driven by a spectrally 
negative stable process, provided a good fit to quoted option prices across maturities. 

By way of application, we employ the sextuple law to obtain semi-analytical expressions for 
the expectations of a number of path-functionals of the drawdown process of an exponential Levy 
process, which provide a description of different aspects of the riskiness of the model: (i) the 
probability that, on a given time horizon, a new minimum is attained (e.g. by a jump) at the first 
moment of a drawdown of a given size; (ii) the expected size of the drawdown process at the first 
moment that a drawdown of a given size occurs, given that this happens before a finite time- horizon 
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and (iii) the probability that, on a finite time-horizon, a drawdown of a given size occurs before 
a raUy of a given size. These and related explicit expressions can form the basis for the analysis 
of optimal investment problems involving the drawdown process — in the interest of brevity, such 
investigations are left for future research. 

In the literature two approaches have been successfully adopted to identify the distributions 
of path-functionals of the drawdown process: the martingale approach, exemplified in e.g. |10^ 
which rests on the identification of certain martingales, and the excursion-theoretic approach, 
employed in a.o. [31 [U [TTl [12l [231 [31] , which is based on the theorem by Ito stating that the process 
of excursions of Y away from zero forms a Poisson point process. In this paper we will follow the 
latter approach. 

The remainder of the paper is organized as follows. In Section [2] preliminaries are reviewed and the 
notation is set. In Sections 13.11 and 13.21 the law of the sextuple is derived for a spectrally negative 
and spectrally positive Levy process, respectively. In Section [Jj the sextuple law is employed to 
derive analytically explicit identities for a number of characteristics of drawdowns and rallies in an 
exponential Levy model. The proofs of the main theorems are contained in Section [5l 

2. Preliminaries 

In this section we set the notation. For the background on the fluctuation theory of spectrally 
negative Levy processes refer to Bertoin [6l Chapter VII] and Kyprianou |22j . 

Let X = {Xt)t>o be a Levy process defined on ((7, J^, (J^()t>o, P), a filtered probability space 
which satisfies the usual conditions. Assume throughout the paper that X is spectrally negative. 
Under P we have Xq = P-a.s. For any x G M we denote by P^ the law of the Levy process started 
from X, i.e. the law of the process x + X under P. To avoid the case of trivial reflected processes we 
exclude X that has monotone paths, i.e. X is assumed to be neither the negative of a subordinator 
nor a deterministic drift upwards. Since the jumps of X are all non-positive, the moment generating 
function E[e^'''"'] exists for all > and is given by il^{9) = t~^ logE[e^"'''*] for some function il^{0). 
The function ip which is well defined at least on the positive half-axis where it is strictly convex with 
the property that limg-j^oo V'(^) = +cc. Let $(0) be the largest root of ^^{O) = 0. On [<I>(0),oo) the 
function ip is strictly increasing and we denote its right-inverse function by ^ : [0, oo) — ?• [<I>(0),oo). 

For q > 0, there exists a continuous increasing function W^"^^ : [0, oo) — )■ [0, oo), called the q-scale 
function, with Laplace transform 

/•oo 

(2.1) / e-^''W'^'i\x)dx = {i;{9) - q)~\ 6 > ^{q). 
Jo 

The function W^''^ is extended to x € (— oo,0) by W^'^\x) = 0. A related g-scale function Z^"^^ is 
defined by 

(2.2) Z^^Xx) = l+q [ W''''\z)dz, for x G M. 

^0 
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The g-scale function W^'^^ is left- and right-differentiable on (0, oo) and we denote the right- and 
left-derivative of W^'^'^ by W^^' and W^^' respectively for any q > 0. Furthermore if a Gaussian 
component is present then for any q >0 we have W^'^^ S ^^(0, oo) (see |13j). In this case we denote 
by W^'^'^' and W^"^^" the first and the second derivative of the scale function W^'^^ respectively. 

We next briefly review two-sided exit results which will be employed in the sequel. Let u < v 
and X E [u, v] and define the first-passage times T~ and T+ as 

T~ = inf{t >0:Xt<u} and r+ = inf{t >0:Xt>v}. 

Let Tu^v = A be the first time X started at x enters the set M \ [u,v\. The two sided 
first-passage results [7] read 



(2.3) 
(2.4) 



-qTu,v J 



VF(9)( 



X — u) 



W(i){v-u)' 
Z^'i\x - u) - Z^'i\ 



X G [u, v\ 

W^i^{x-u) 



V — u) 



WM{v-u)' 



X G [u, v\. 



A Levy process started at zero creeps downwards (resp. upwards) over a level x < (resp. x > 0) 
if, with positive probability, the first time that it enters the set (— oo,x) (resp. (x,oo)) this does 
not happen by a jump. It is an immediate consequence of the Wiener-Hopf factorization that a 
Levy process creeps both upwards and downwards if and only if its Gaussian component is not zero. 
Therefore a spectrally negative Levy process with non-monotone trajectories creeps downwards if 
and only if it has a positive Gaussian coefficient (e.g. [H p. 175]). In this case the results of 
Millar [28] imply that 



(2.5) E, 



— (w^''^'{x-u)-^iq)W^i\ X — u)) for any u < x, q > 0. 



Both sides of the equality in (j2.5p are understood to be equal to zero if o" = 0. The formula for 
the probability that X leaves the interval [u,v] by hitting u follows from ()2.3p . (j2.5p and the strong 
Markov property and is given by the following expression: 



(2.6) 



E,: 



a 

T 



X — u) 



W^'i'^'{v-u) 



1^(9) (W - U 

for X G [n, v\. Again the expression is understood to be equal to if o" = 0. 



X — u] 



3. Reflected spectrally negative Levy processes 

Let Xt = supo<u<t {Xu] and = info<M<t {^il} and define the reflected processes Y = {Yt)t>Q 
and Y = {Yt)t>Q by 



(3.1) 



Yt = Xt-Xt and Yt = Xt - X^. 
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The focus of this paper are the first-passage times over a level a > of the reflected processes Y 
and Y, 

(3.2) Ta ■■= mf{t >0:Yt> a} and 9a ■= mf{t >0:Yt> a}. 

It is well known that the stopping times and Ta are finite P-a.s. Furthermore, in [3^, Theorem 1] 
the Laplace transform of was identified as 



(3.3) E[e-''"'''] = Z^i^c 



The last times before time t that X visits its running supremum and infimum are denoted by 
and respectively, where 

(3.4) Gt = sup {s <t : Xs or Xs- = Xs} and = sup{s < f : or = X J. 

Note that Gt (resp. can be viewed as the last time before time t that the reflected process Y 
(resp. Y) is equal to 0. 

In Section [3.11 we characterise the joint law of the following sextuple of random variables: 

Ta the first-passage time over a level a of the reflected process Y, 

Gra the last time that X was at its running supremum prior to the first-passage time Ta, 

Xra the supremum of X at the first-passage time Tq, 

X_^^ the infimum of X at the first-passage time r^, 

Yra- the position of the reflected process just before it crosses the level a, 

Y-r^ — a the overshoot of the reflected process Y over the level a. 

In Section 13.21 we give the joint law of the following quadruple of random variables 

Ta the first-passage time over a level a of the refiected process y, 

G-^ the last time that X was at its running infimum prior to the first-passage time Tq, 

the supremum of X at the first-passage time Tq, 
Xy^ the infimum of X at the first-passage time Tq. 

Note that in this case, since X is assumed to be spectrally negative, the refiected process Y can 
only jump down. Since Y is right-continuous with left limits at the first-passage time Ta we have 

„ Qj ct.S. 

'a 'a 

3.1. The sextuple law. We now give the law of the sextuple (r^, Gt-„, X^-^, X^^, 5^^^-) — a). 
Define for any a > and p, > the map Fp^q^a ■ IR+ — ^ 1^+ by 

(3.5) Fp^<y,a(?/) = A(a,g)exp(-yA(a,p)) , y G K+, 

where X{a,q) is the ratio of the derivative of the g-scale function and the g-scale function at a, 

(3.6) A(a,,) = ^5$^. 
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Let A be the Levy measure and > the Gaussian coefficient of X. Denote by R^\dy) = 
E [Jq^" e~''*I{yjg(jy}dt] the g-resolvent measure of Y kihed upon first exit from [0, a], which may be 



expressed in terms of the g-scale function VF^'?) by ([301 Theorem 1]) 



(3.7) 



R^^\dy) = \xia, W^'^\dy) - W^'i\y)dy 



y G [0, a] 



Furthermore, let A(«)(a) denote the expression that is given in terms of the derivatives of W^"^^ by 



(3.8) 



A('')(c 



a 

~2 L 



W^'^\a)-\{a,qY^W^'^"{L 



where the expression on the right-hand side of the equality in (j3.8p is taken to be zero if o" = 0. As 
we will see below, the Laplace transform of on the event that Y creeps over the level a is equal 
to AW(a)=E[e-''-'^/{y^^=,|]. 

Theorem 1. Let Xq = x G M and a > and define events 

= {X^, > u.'Xra e du, y^„_ E dy, Yr^ - a £ dh} and = {X^^ > u,Xra e dv, Y^, = a} , 
where u,v,y and h satisfy 

(3.9) u < X, y £ [0, a], v > xy (u + a) and h £ {0,v — u — a]. 

Then for any q,r >0 the following identities hold true: 



E,: 



-qra-rGra 



(3.10) 



(3.11) 



E. 



w^'i+''\{x - u) ^ . 

Ty('?+'^)(a) 
W^i+''){a) 



■ Fq+r,q,a {v - {x V {u + a))) R^^^ (dy) K{y-a-dh) 



q+r,q 



,a iv-{xV{u + a)))A^i\a) 



where /|.} denotes the indicator of the set {•} and cf\d = min{c, d}, c\/ d = max{c, d} for c, d G R. 

Remarks, (i) The formulas in (|3.10p and (j3.1ip determine the law of the sextuple 
[Ta, Gr^, Xt^, 2^^^,Yt-^-,Yt-^ — a). Indeed, note that if the parameters u,v,y and h do not sat- 
isfy the restrictions in (|3.9p . then we have P^,. [Aq] = [Ac] = 0. In particular, if we take u = — oo, 
then ()3.9p places a restriction neither on the size h £ (0, oo) of the overshoot of the reflected process 
nor on the level u > x of the supremum of X attained at Ta- This is also transparent from Figured! 



(ii) Note that on the event Aq we have X^^ = v — a — h Pa;-a.s. and hence X^^ = {v — a — h) AX 
Define the event 

Ao = {X^,- > u,Xr, £ dv, Yr,_ £ dy, Y^^-a£ dh] . 
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Then the following implications hold: 

h£ {0,v -u - a] =^ A~ = Ao ¥^-&.s. 

h > V — u — a =^ X^_^ = V — a — h Fx-a..s. on A~ and P^,. [Aq] = 0. 



Furthermore, it is clear from step (2) of the proof of Theorem [T] (presented in Section [5|) that 
Ej; e"'''^""^*^^''/^- is given by the formula in (I3.10p . which in this case holds without a restriction 
on the size of the overshoot h from (13. 9p . 

(iii) Note that the reflected process Y creeps over the level a (i.e. P[1V„ = a] > 0) if and only if X 
creeps downwards which is the case precisely if its Gaussian coefficient cj^ is strictly positive. Hence 
the expression A^'^\a) given in (|3.8p is well-defined, as W^'^'^ is twice differentiable at a > if cr^ is 
strictly positive. 

(iv) Schematic representations of a typical path of the process X in the events Aq and A^ of 
Theorem [1] are given in the bottom and top pictures in Figure [TJ Note that Ac contains the paths 
of the refiected process Y when it creeps over the level a and that l^a- = « P-a.s. on Ac- Similarly 
Aq consists of the paths of Y which first enter (a, oo) by a jump. 

Explicit expressions, in terms of the 0-scale function, for the marginal distributions of X-^a and 
X^ under P^, can be obtained from Theorem [Tl 

'a 

Corollary 2. For any a > and x S M, the following hold true: 

(3.12) Fx{x-Xra~<z) = W{zAa)/W{a), z G M+, 

(3.13) Px(Xra -x>z) = exp {-zW!^{a)/W{a)) , z G M_|_, 

where W = W^'^^ denotes the 0-scale function. 

Proof of Coro//ary [H Setting q and r equal to and letting u = — oo in ()3.10p - ()3.1ip . and integrating 
over y £ [0,a], h £ (0,oo) and v G (x,oo) the right-hand sides of the these expressions we obtain 
that c(a) = 1 where 

poo pa 

(3.14) c(a) := / / i?i°)(dy)A(y - a - dh) + AW(a). 

Jo Jo 

Fixing x,a,Uo,Vo, such that a > 0, uo + a < vq and x £ [uq, vq], and again integrating the right-hand 
side of (|3.10p over y £ [0,a], h £ (0,oo) and v £ {vo,oo) and the right-hand side of Eqns. (j3.1ip 
over V > Vq we find 

P.(^..- > no,X.„ > .0) = c(a) • ^^"^^^w";/''"^ e-(-o-v(«o+a))AM)^ ^ ^ ^. 

We obtain ()3.12p and ()3.13p from this display by setting vq = xV{uo+a) and uq = — oo, respectively, 
and recalling that c(a) = 1. □ 
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Figure 1 . The top fi gure contains a schematic representation of a typical path of the process X 
in the event Ao = {]Lt — '^i ^ dVjYr^- £ dy,Yr^ — a £ d/i}. This figure illustrates the idea 
behind the proof of Theorem[T] the depicted path satisfies Tu,u+a = T^+a ^^'^ ~ Gra- It is clear 
from the figure that the trajectory can be decomposed into the following three parts, each of which 
is analysed separately in the proof of Theorem[T] (1) the segment over the interval ^,Tu,u+a\, (2) 
the jump at time Ta and (3) the segment that straddles Gt„ over the time interval (T^.u+a, Ta). The 
bottom figure depicts a schematic representation of a typical path of the process X in the event 
A, = {X^^ > u,Xr^ e dv, n„ = a}. 
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Remarks, (i) If X is a Brownian motion with drift started at zero, i.e. Xt = fit + crWt, where 
{Wt)t>o is a standard Brownian motion and /i G M, > 0, then the scale function takes the form 

(3.15) W{x) = -h-e'^""] if 

/i L J 

(3.16) W(x) = 4tX if /i = 0. 

Since X has no jumps, we have an almost sure equality Xr^ = a + X^^. Therefore the random 
variable a + X^-^ is by p.l3p exponentially distributed with parameter (2/i/cr^)/(e°^'^/'^^ — 1) (resp. 
1/a) if /i ^ (resp. /i = 0). This fact, first observed by Lehoczky [23], is generalised by the formula 
in p.l3p to the class of spectrally negative Levy processes. 

(ii) Let X be a spectrally negative a-stable process with cumulant generating function -0(0) = [a6)°' 
for > 0, where a G (1) 2], a > 0. Identity (|2.ip and the definition of the Gamma function T imply 
that the scale function of X takes the form 

(3.17) W{x) = x>0. 



Corollary [2] implies that, under P^,., the random variables a ^ yx — X^ j and a ^ (X^-^ — x) follow 
Beta(a — 1, 1) and Exp (^^) distributions with respective probability density functions 



^.{x - X,^_ e dz) = /[o^,](z)^_ (^-j dz, 

P,(X.„ - .X G dz) = /(o^^)(z)^e-^("-i)/'^dz. 

Note that these distributions do not depend on the value of cj. In the case of Brownian motion 
started at x (i.e. for a = 2), the random variables (x — X^J and o"^ {Xr^ - x) follow [7(0, 1) 
and Exp(l) distributions respectively. 

3.2. The quadruple law. In this section we characterise the law of the quadruple 
(tq, Xf^), where G^, defined in (j3.4p . is the last time the process X was at the infimum 

prior to time t. Before stating the result, we recall that q i— )• W^'^\x) has a holomorphic extension 
to C for every x > and that (x, g) i— >• Z^'^\x) has a continuous extension to [0,oo) x M, which is 
holomorphic in q for every x G M (see [30], Lemma 2). Furthermore for any u > let P" be the 
exponentially tilted probability measure defined via the Esscher transform (see [H Ch. XIII]). Then 
the scale functions x i— )• Wu ^^^'^\x) and x i— )■ Zu '^^"^^(x) satisfy the following identities for all 
g G C,x G M: 

(3.18) Wi'^-'''("»(x) = e-"^Vr('?)(x), 

(3.19) ^^^-^(""(x) = l + {q-ij{u)) r e-'''W^'^\z)dz. 



The identity in (j3.18p . for q > ip{u), follows by taking Laplace transforms on both sides and 
applying ()2.ip . and hence by analyticity for all g' G C. The identity in p.l9p follows from ()3.18p 
and the definition in (12.21). 
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V — a 




Figure 2. A schematic representation of a path which satisfies Tv-a,v = T~_^ and is therefore 
in the event {Xf^ < v}. At the moment of local time t the path in the figure has an excursion et 
away from the infimum. The process is the inverse local time at zero for the Markov process 
Y. 



We now state the main result of this section. 
Theorem 3. Let q, r,u,v > 0. Denote p = q + r — ip{u). Then the following identity holds: 



(3.20) 



1^(9+0 (a) 
Wi'i){a) 



Remarks, (i) li a < v, then the definition of imphes P \_X^^ < v\ =1. Therefore Theorem [3l 
together with equations (|3.18p and (|3.19|) (recall that W'^'^\y) = for all (7 > and y < 0), implies 
the identities 



(3.21) 
(3.22) 



E 



-qTa+uX_~ 



Efe- 



-QTa-rG- 



^(9+0 (a) 



g, li > 0, 
g, r > 0. 



' Ty(9)(a)Z(9+0(a)' 

The special case of formulae p.2ip and (I3.22p for u = and r = respectively (i.e. the Laplace 
transform of To) is well-known (see [30j . Proposition 2). Note that in (j3.22p the random variable 
under the expectation does not depend on the starting point of the process X. 
(ii) Since X is spectrally negative, we have 1%^ = a almost surely. Hence formula (|3.20p also yields 
the joint law of the quadruple (ra,G:^^Xf^,Xf^). 
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(iii) Under the law P^, the cdf of the random variable Xi^^ — x is given by 
(3.23) p,(X,.-.<.) = l-M^Z_£lM, ,eR. 

Note that the distribution of X^^ — x under P^. is absolutely continuous with respect to the Lebesgue 
measure if and only if 1^(0) = 0. If VF(0) > 0, this distribution has an atom only at a, which is 
of size W{Q)/W{a). Furthmore, the Laplace transform of x — X_^^ is given in terms of the scale 
function of X by the formula 



(3.24) E^. \(r<''-^-a)] = e-"7 



l-i^iu) / e-"^Ty(z)dz 







(iv) In the special case when X is a spectrally negative a-stable process with cumulant gen- 
erating function 'ip{6) = {(jO)", where a £ (li^], a > 0, the scale function takes the form 
W{x) = x"~V(c^"r(Q)), for X > 0, by (IXTTll . Formula in ([3:23]) implies that under P^;, for any 
X G M, the random variable {^fa ~ ^) follows a Beta(l,a — 1) distribution with probability 
distribution 

F^{X^^-xedz) = /[o,^](z)^— dz. 
The Laplace transform of the random variable (x — X^^) is, by ()3.24p . equal to 

(3.25) E Je~"'^"H^-2^?Jl = e-"-lM^, where T{z,y) = T e-'s'~^ds 
L J T{a,u) Jy 

is the incomplete gamma function. Note that the laws of {^Ta ~ ^) ^^'^ "^^^ ~ —fa) d° 
depend on the value of a. In particular, in the case of Brownian motion started at x (i.e. a = 2), 
the random variables {^?a ~ ^) ^^"^ '^"^ ~ ^9a) fo^ow f7(0, 1) and Exp(l) distributions 
respectively. 

4. Application: drawdowns and rallies in exponential Levy models 

Let the stochastic process S = {St)t>o model the price of a stock or a foreign exchange rate. The 
(absolute) drawdown and rally processes of S are defined by the difference 5 — 5 of its running 
supremum St = supo<„<( {Sa} and the current value, and the difference S — S_oi the current value 
and the running infimum S_t = info<M<t {5'u}. Their relative counterparts, the relative drawdown 
process and relative drawup or relative rally processes of S are given by S/S and S/S_. Let Dq, 
be the first time the price process S drops below its running supremum by at least (100a)% with 
a G (0, 1), and denote by C/^ be the first time the stock price S rallies above its running infimum by 
at least (100/3)% with /3 > 0. In some trading strategies buy- and sell-signals for S are generated on 
the basis of its relative drawdown or drawup processes. For example, a commonly used strategy is 
to buy at the epoch C/g, and to sell at the epoch Da (see e.g. Levich &; Thomas [25] for an analysis 
of such a strategy in an foreign exchange setting) . In this section we will consider the following four 
criteria that provide descriptions of different aspects of the risk associated to investing in S: 
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(i) The probability that S attains a new (ah-time) minimum at the moment Da of a first 
relative drawdown of size (100q;)%, on the event that Da is before T, where T > is a given 
time-horizon. 

(ii) The probability that S attains a new maximum at the moment {/^ of a first relative drawup 
of size (100/3)%, on the event that Up is before T. 

(iii) The expected (absolute) drawdown of S at the epoch D^, on the event that Da is before T. 

(iv) The probability that a relative drawdown of size (100a)% occurs before a relative drawup 
of size (100;5)%, on the event that Da is before T. 

Assume henceforth that the price process is modelled as 5 = {St)t>o where 

St = Soexp{Xt), t>0, 

with 5*0 > and X = {Xt)t>o a completely asymmetric Levy process. Then, using the notation of 
Section[3l we have S = , S_ = e^, and therefore S/S = e^, S/S_ = e^, and find 

(4.1) Da = inf {t>0 : St/St > 1/(1 -a)} =Ta, where a = -log(l-a), 

(4.2) U/s = mf{t>0:St/St>l + (3} = Tb, where 6 = log(l + /3). 

Furthermore, the quantities in (i)-(iv) can be expressed in the notation of Section [3] as follows: 

F[SD^=So^,Da<T] = ¥[Ta = G,^,Ta<T], 

nSu, = Su,,Uis <T] = F[Tb = Gr„n<T], 
n(SD^-SDjI{D^<T}] = 5oE[(e^--e^-)/|,„<r}], 

F[Da<Up,Da<T] = F[Ta<Tb,Ta<T]. 

In the following sections these quantities are explicitly expressed in terms of scale functions, em- 
ploying the sextuple and quadruple laws that were derived in Section [3l 

4.1. The probability of attaining a new minimum at the first moment of a drawdown. 

In the next result the distributions of the two random vectors (Ta,l^„) and (ra,!^^) are explicitly 
identified in terms of scale functions. Note that the process X is equal to its running minimum at 
the first moment of drawdown of size a if and only if is equal to zero. Similarly, X is equal to 
its running minimum at the first moment of a drawup of size b if and only if Y^^ is equal to zero. 

Corollary 4. Let a > 0. (i) For b,q,6 > the following identities hold true: 



(4.3) 
(4.4) 



E 



E 



{Yra>b} 



-bX{a,q) 



(g_^(g))e-^ 



W('i){a) 



WM(a) / Jo 



M, 



q,ai 



VF(9)(z)dz, 
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where X{a,q) = w!^^' {a)/W'^''\a) and Mg^a ■= E[e~^('''^)(^-^--'')] is given by 



(4.5) M,,a = z[-^^^'^\a)-W. 



with A = A(a, q). 



(a) The Laplace transforms of t ^ P(ra < t,Ta = G^^) and 1 1— >■ P(ra < t,Ta = G^^) are given by 

/•oo ^ I 

(4.6) 



^"e-^*P[ra<t,r, = G,Jdt = ^|l-A(a,g)^ 



1-X(a,q) /° dy\M„a, 



(4.7) 



— if r w^'^Hv) 1 



Remarks, (i) From the formulas in (j4.6p and (j4.7p it follows that the non-negative random variables 
Yr^ and have atoms at zero of the sizes 

w'^{a) r_ 

10 



(4.8) 
(4.9) 



^r„ = 



1- 1 



w' (n) 



Furthermore, the formula in (|4.3p . implies that conditional on the event > o|) follows an 
exponential distribution with mean W{a)/W'^{a). 

(ii) In the case X is Brownian motion with non-zero drift fi and Gaussian coefficient <t^, the scale 
function is given by (j3.15p . Hence the atoms of and at zero are by (j4.8p and ()4.9p of the 
sizes 



i^r. =0 



e«2M/<x2 _ 1 _ a2/i/cj2 



Furthermore, ^-^id Y^^ conditioned to be stirctly positive follow exponential distributions: 



y^jy^„ >o~Exp( 



> ~ Exp 



2/z/cj2 



If /i = 0, formula (|4.8p implies 



Y^^ = 



1/2. In this case formula ()4.3p implies that 



P G d6|y^„ > = (e-''/"/a)d6 for 6 > 0. These results coincide with the findings of [E 
Prop. 2.2, Cor. 2.6 . 

(i) For later reference we note that the formula in (j4.4p can equivalently be expressed as 
(4.10) 



E 



e-^'Wl'>'{a + z) dz r W^'i\z) 
^ ~ /o°°e-^-iy('?)(a + z)dz io 



for any g > and > ^{q). 



Proof of Corollary ^^The identity l^a = — ~ ^Lra implies that, for any /i > 0, on the event 
Yra <b,Yra-a£ d/i} = {Xr^ < a + b + h + X^^,Yr^ - a £ dh} 
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the inequality < a + b + h holds Po-a-s. The formula in ()3.10p of Theorem [T] and the definition 
of -Fq,g,a given in (j3.5p imply 



(4.11) 



E 



-qTa 



{Yr^<b,Yra-a€dh} 



h{dh)+l2{dh), 



where Ii and I2 are measures on (M+,;B(M+)) given by 



Ii{dh) := E 



{Yra <b,YTa ~aedh,Xra <b+h} 



-QTa J 



{Yra <b,yTa -aedh,X e{b+h,a+b+h] } 



Wi^Haj 



dy. 



In these expressions \{a,q) is given in dSSl) and 4a, i G {0,1}, are measures supported on M+ and 
defined by the following formula 



{dh) = e-^^^'--'"^ / (dy) A (y-a-dh), 



(4.12) 



where is the resolvent measure given in ()3.7p and A is the Levy measure of X. An analogous 
argument, based on the formula in (I3.11|) of Theorem [U yields the formula 



(4.13) 



E 



{Yr^<b,Yra=a} 



1 - e-^^('^'^) + A(a, g)e-^^("'^) ^^^^ dy 
^ ' Jo W(i){a) 



A(^)(a), 



where A^^\a) is given in (j3.8p . Observe that, on account of the compensation formula applied to 
{AXt)t>o, the following relations hold true: 



(4.14) 
(4.15) 



^[6-"^"] = i^°a(0,oo) + A('?)(a), 
E[e-^(^^»-)] = (O,oo) + A(^)(a), X = X{a,q), 



where the expressions for E[e ^"^"J and E[e in terms of scale functions are given in (j3.3p 

and (|i3]) (see [31 Thm. 1]). By integrating (|4.1ip over h £ (0,oo), adding (|4.13p and using (|4.14p . 
(j4.15p and (|3.3p we obtain the formula in (|4.3p . 
It follows from Theorem [3] that the identity 



(4.16) 



E 



p-QTa+uX^^j 



z('^"'/'(«))(a-^;) 



a — V 



(a) 



holds for q,u,v > 0. The Laplace transform in v of the identity in ()4.16p evaluated at n = ^ implies 



the equality in (j4.4p . Since E 



limg^oo E [e" obtain the formula in 



(j4.7p by noting that {1^^ = 0} = = G^^} and applying the identity 



lim 



9^00 0/°°e-^^H^(5)(a + z)dz VF(5)(c 



which follows from well-known properties of the scale functions. 



□ 
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4.2. The expected drawdown S — S at Da- 

Corollary 5. The Laplace transform ofT^ E[(S/)^ — SDa)^{Da<T}\ ^-^ given by 



e-'^''E[{SD^-SDjl{D^<T}]dT 



(4.17) 



SqX 
A - 1 



hg{x)W^''\x)dx 



So 
A- 1 



hg{x)W'''^\dx) 



where a = -log(l - a), A = X{a,q) = w!f'' {a)/W^'i\a) and hg{x) = 1 - (1 - g-V(l))e"'^. 

Proof. Recall that a = — log(l — a). In the notation of Section |3l the Laplace transform on the 
left-hand side of equation (14.17P is equal to 

— E[e-«^°|e^-- - e^"-|l = ^E[e~«^"-+^"''|e-'?("'"-^"'")(l - e"'^"")!!. 
q q 

Direct integration of the expressions in Theorem [T] with u = —oo and x = shows that, for q > 

(4.18) E[e"''"''"{e^"- -e^"-}] 

= E[e-«^-«+^^"]{E[e"^("'''"^"'")] -E[e-«(^''-^--)-'^^-]}, 



with 

^^Q-<l{-ra-Gra)-iYTa-a)j 

(4.19) 

(4.20) E[e-«(""°-^-«)] 

(4.21) E[e-''^"-+^"-] 



A(Q,g) 
A(a,0) 



Z 



a) - W, 



'(a) 



1+ 



(a) + 



(9-^(1))/ 



X{a,q) 



Z^^^a) -qX{a,q)-^W^'>\a) 



A(a,0) 
A(a, 0) X{a,q) 



X{a,q) X{a,q)-V 

where the functions zj"^ ^^^^\a) and W^'^ ^^^^^ are defined in equations p.l9p and p.lSp . In- 
serting the expressions in equations ()4.19p and ()4.20p into equation (I4.18P yields the expression in 
equation (11171) . □ 

4.3. Probability of a large drawdown preceding or following a small rally. In this section 
we identify explicit expressions for the probabilities of the events {Da < Ujs} and {Da > f/^g} in the 
case when the size 1/ ( 1 — a) of the relative drawdown is larger or equal to the size 1 +/3 of the relative 
drawup. Note that the case 1/(1 — a) > 1 + /3 is the case of most interest in practice, since the 
event {Da < f^/?} corresponds to the risk held by an investor of losing more in an investment than 
they would gain if the security rallied. We refer to Hadjiliadis & Vecer [19] for further background 
on the use of this probability as risk-measure. 

The results in this section are a consequence of Corollary H] and the following key identity: 



16 



ALEKSANDAR MIJATOVIC AND MARTIJN R. PISTORIUS 



Lemma 6. Pick any numbers a > b > 0. Then the following equality holds almost surely: 
(4.22) {Ta<n} = {a-b<Y^^}. 

Remark. In the case b > a and the Levy measure of X is non-zero, the equahty analogous to (|4.22p 
does not hold. The event {r^ < Tf,} is different from {b — a > Yr^} as it depends in an essential way 
on the supremum Yr^ of the reflected process Y at r^. The identification of the law of the random 
variable Y^-^ is beyond the scope of the current paper. 

Proof of Lemma\^ Let Y and Y denote the processes given by the running suprema of the reflected 
processes Y and Y respectively. In Hadjiliadis & Vecer [T2] the following model- free equalities were 
established: 

'ovfy-y 



(4.23) Y + Y = Y yY = Y + 

Since Y^^ = Y^^ = b a.s., we find {tq < = {a < Y^^} = |a — 6 < Y^^ — Y^^^, and the lemma 
follows in view of (j4.23p . □ 
By combining Lemma[6]with Corollary[5]we obtain a semi-analytical expression for the probability 
of interest: 

Corollary 7. Let a > b > and t > 0. Then we have P < A t] = P [Y^^ < a - b,Th < t] , 
where the joint Laplace transform of (u, t) i— )• P [l^^^ <u,%< t] is given by the formula 



poc poo 

/ / e~'^-i'F[Y^^<u,?b<t]du 
Jo Jo 



(4.24) ^ 



\^(m^^ o__\ rV(^)(,)d, 



for any 6,q > 0. Moreover for any a > we have 

W[{a) 



W a) r"- 

1^°<^J = WW I 



Remarks, (i) It is a consequence of the representation of the scale function Zg '^'^^^'^ given in ()3.19p . 
that the Laplace transform in (j4.24p can be expressed as 

e-^"P [Y^^ <u]du = l+ 1^1/ J e~'^W{b + z)dz- 9/W{b)j J W{y) dy. 

(ii) The probability measure in Corollary [7] does not depend on the starting point of the process 
X, since the reflected processes Y,Y, and hence the stopping times Ta,%, are independent of Xq. 
Therefore equality in (j4.24p remains valid if P[-] is replaced by Px-['] for any starting point x € R. 

(iii) Recall from p.lSp that in the case X is a Brownian motion with drift G M \ {0} and the 
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Gaussian coefficient > 0, the scale function takes the form W{x) = ^1 — e~^^^'/°"^^ /fj, for x > 0. 
The formula in Corollary [7] yields 

e-W^- - 1 + a2^,y 

which coincides with the formula in (2.1) of |I9]. In particular if = then Corollary [7] yields 
IP [7"a < ^a] = 1/2, which follows of course also directly by symmetry. 

4.4. Example: Carr & Wu model. Carr &: Wu [9] documents the persistence of the implied 
volatility skew across maturities in the exchange traded options on the S&P 500 index. The evidence 
presented in [9] suggests that the left tail of the risk-neutral distribution of returns of the index 
remains "fat" as the maturity increases. In order to capture this phenomenon Carr & Wu [9] model 
the risk-neutral evolution of the index level by the following stochastic differential equation: 

(4.25) dSt/5t_ = (r-g)dt + SdL"'~\ qG(1,2), S>0, 

where r and q denote the continuously compounded risk free rate and dividend yield. The driver 
L"'~^ is assumed to be an a-stable Levy martingale with maximal negative skewness, i.e. L"'~^ 
is a spectrally negative a-stable process. In this model the decay of the left tail of the log return 
(i.e. of the quantity P [log (St /-St) < s] for any fixed < t < T) is asymptotically equal to 
as s tends to — oo, while the right tail decays exponentially, so that the model (j4.25p captures the 
observed phenomenon of the persistence of the implied volatility skew across maturities. 

In order to understand the risk-neutral probabilities of drawdowns and rallies in this model for 
the S&P 500 index, note that the process S can be expressed as S = exp(X), where X is a spectrally 
negative a-stable process with drift. The cumulant generating function of X takes the form 

(4.26) tpie) = fiO + [aO f for 6 > 0, where a G (1,2), 

for some G M, o" > that depend on r, q, S and a. Note that in the limit case of a = 2, the process 
X corresponds to Brownian motion with drift and variance 2(7^ at time one, and in this case the 
model ()4.25p reduces to the Black-Scholes model. 

As can be seen from Proposition [71 the probability of a drawdown preceding a rally in the 
model (j4.25p is given in terms of scale functions of the process X, which is expressed in terms of 
the Mittag-Leffier function 



oo 



,,n 



Note that -E/s,^ is an entire function if /3, 7 are strictly positive and that it generalises the exponential 
Ei,i{y) = e^ and the hyperbolic cosine i?2,i(y^) = cosh(y) functions. Furrer jl8j identified the 0-scale 
function 14^ of X as 



1 



(4.27) W{x) = - 1-E„_i,i -^x"-M , rE>0,aG (1,2]. 



a"- 
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As a consequence, the g-scale function W^'^'^ of X, which is related to W by the weh-known expression 
W^'^^x) = Y.k>o^^^*^^^^H^) that is given in terms of the A:-th convolution W*'', k e 'N, of W 
with itself, admits the following series representation: 



(4.28) 



W'^''\x) = - 



l-E, 



0-1,1,1 



where, for any nonnegative /?, 7 and S, the function Ep^^^^ : 

by 

00 00 



n/3 + k5 + 7) 



x,q>0,ae (1,2], 

M is the function that is defined 

y,zeR, 



n=0 k=0 

where the coefficients {sk,n, n,k £ NU {0}) are given by 



Sfc,n+1 



n 
k 



/c,n G NU{0}, 



so,o = 1- 



Note that Sk n 



k\{n-l~ky. 



if = 0, 1, . . . , n — 1, and n is equal to zero otherwise. 



Remarks, (i) In the limit as a approaches 2, the formula in (14.270 tends to the scale function of the 
Brownian motion with volatility (J\/2, which is given in ()3.15p . 

(ii) It follows from the formula in (j4.27p that if the modulus of the drift /x tends to zero, then in 
the limit we obtain the scale function, of the spectrally negative a-stable process, which is given 
in (131711 . 

(iii) It is assumed in (j4.27p that a > 1 and hence the series Ea-i^i converges uniformly on compact 
subsets. Therefore the derivative and the integral of the scale function can be obtained by performing 
these operations under the summation: 

a — 1 



(4.29) 
(4.30) 



W'{x) 
W{z) dz 



-X ^a~l.l 



a 



l-E, 



a-1,2 



X > 0. 
a > 0. 



We now present an analytically tractable expression for the probability of a large drawdown 
occurring before a small rally. Let denote the inverse Laplace transform. 

Corollary 8. Let a = — log(l — x), x £ (0, 1), and b = log(l + y), y > 0, and a > b. Then the 
probability of a drawdown of size (100x)%, in the Carr & Wu model (I4.25P occurs before a rally of 
size (100y)% is given by P [Ta < %] = C-^{F){a - b) where 



1 



(4.31) F{e) = - 



/xe 



/ E 



\n=\ 



afla— 1 



r(n(a - 1) + 1,66') 
r(n(a - 1) + 1) 



+ 1/VF(6) 



W{z)dz, 



whereT{n{a — l) + l,b6) denotes the incomplete gamma function (see (|3.25p ) and the scale function 
W and its anti- derivative are explicitly given in (|4.27p and (j4.30p . 
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Remark. Standard Laplace inversion algorithms, such as the one given in Abate and Whitt [1], would 
typically evaluate the function in (j4.3ip about twenty times in order to compute the probability 
P [Y^^ < a — b\. Since the series in expression (j4.3ip is dominated by a geometric series, the error 
at each evaluation of the Laplace transform can be controlled easily. 

In the case when the relative drawdown level 1/(1 — x) is strictly larger than the size 1 + y of the 
relative drawup (i.e. rally) of the index S, the probability of the market crash P [ra < where 
o = — log(l — x) and b = log(l + y), is by Proposition ([7]) equal to P [l^b < a — b\. By Proposition [7] 
the Laplace transform of the function u i— )• P [l^^^ < u] is given by the formula in (j4.3ip . 

The previous results yield a closed-form expression for the probability of a drawdown occurring 
before a rally in the case these are of the same size in the log scale (cf. (j4.ip and (j4.2p ): 



Corollary 9. The probability of a drawdown of size (100x)%, x G (0,1), in the Carr & Wu 
model (j4.25p occurs before a rally of size (100x/(l — x))Vo is given by 

(4.32) P [Ta < 9a] = {a- l)AE' [-A] ^} ~ f""^''/"!, , with A = -^a^~\ 



:i-i?,_i,i [-A]Y 



where a = — log(l — x). 



Note that in the limiting case a = 2, formula (|4.32p yields the well-known expression for Brow- 
nian motion with drift ji and volatility a^/2. Formula ()4.32p follows from Proposition [7] and the 
formulas (g^ZD, and (O0]l . 



5. Proofs 

5.1. Proof of Theorem [1]. We now give a proof of Theorem [1] based on Ito's excursion theory. 

Start by noting that since both sides in ()3.10p and (j3.1ip are continuous in g, it suffices to prove 
the formulae in (j3.10p and ()3.1ip for (7 > 0, which will be assumed without loss of generality. We 
prove the result in three steps. Step (1) deals with the segment of a path in 

Ao = [X^, > u, Xr, e dv, y^„_ G dy, Yr^-ae dh} 

over the time interval [0, r^^u-i-a] by applying the strong Markov property at time Tu,u+a- Step (2) 
extracts the final jump AX^-^ in the expectation of ()3.10p by applying the compensation formula 
for the Poisson point process {/S.Xt)t>o- Step (3) applies the Ito's excursion theory to the segment 
of a path in Aq over the interval {T^^u+a, Ta)- In the case of the event Ac the structure of the proof 
is similar with steps (2) and (3) merged as there is no overshoot at time Tq. 

Step (1) Note that P^^ [{Tu^u+a = T^} n = (see Figured]) and that on Ao we have Tu^u+a < 
Gr, P^-a.s. We can therefore apply the strong Markov property at Tu,u+a 

and identity ()2.3p to the 
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expectation in ()3.10p : 



'qTa-rG-r 



-(q+r)Tu,u+aj 



It _l. -t+ 1-^ 



(5.1) 



VF(''+'')(aA(x-u)) 



E. 



xV(ii+a) 



-qTa-rGr, 



Ia. 



l^('?+'-)(a) 

S'tep ('^J By ()5.ip we may assume that the starting point x of the process X is greater or equal 
to u + a. Since Y^^^ G [0, a] this imphes that 



(5.2) 



X_^^_ > u P^-a.s. and hence {X_^^ > u] = + AX^^ > u} P^-a.s. 



In order to deal with the jump at time we apply the compensation formula (see Section 0.5 
in [6]) to the Poisson point process {/S.Xt)t>o as follows: 



E, 



-qTa-rGra 



Ia. 



E, 



rG 



raj 



E. 



^ g '''^'-^{sup,<i{y4<a,yt>a,Xt6di),yt_edy}^{y-a-AXtGd/i,AXt>M+?/-?^} 
t>0 



I{h&(o,u-v-a)}^ {y-a- dh) E^ 



n-qt-rGtj 



(5.3) 



/{/i6(o,«-i,-a)}A {y -a - d/i) (Ea; X £;) 



^{X^edv,y^Gdj/,r;<r„} 



where the sum runs over all i > such that AXt ^ and, as noted before, q > 0, and rj is an 
exponential random variable defined on a probability space (R+,S(M+),P) such that (P x P)[7] > 
s\J-'] = P{r] > s) = e"'^'^ for s E M+. The second equality follows by (|5.2p . The third equality follows 
from the compensation formula for the Poisson point process of jumps of X and the fact that the 
inequality AXt > u + y — v is equivalent to the restriction on h given in (j3.9p . 

Step (3) The final task is to calculate the expectation in (j5.3p . This requires the Ito excursion 
theory for the refiected process Y . Since X is spectrally negative, the local time at zero of Y can 
be taken to be the running supremum of X, i.e. L = X — x, under P^;. Let (t, et)t>o be the Poisson 
point process of excursions of X from its supremum X, taking values in [0, oo) x (<5 U {5}), where 

S = {ee D{R) : 3C e (0, oo] such that e(C) = if C < oo, e{0) > 0, e{t) > Vt G (0, C)} , 

L'(M) is the Skorokhod space and d is the graveyard state, with the intensity measure dt x n(de), 
where n denotes the ltd excursion measure. For each instance of local time t £ {0,v — x], the inverse 
local time satisfies L^-i = X ^-i —x = t since L^^ = T^_^^. Since the excursions are indexed by 
local time, for any t € {0,v — x] we have ej = {Xj^-i — Xj^-ij^^ : < s < L^^ — L^J:} if L^J: < L^^ 
and = 3 if L^l^ = L'^^ . This implies that the excursions can be indexed by a subset of actual 
time that is given by the left-end points of excursion intervals. For any excursion e £ £■ define 
Ta{£) = inf{s > : e(s) > a} (with the convention inf = oo) and let C,{e) be the life time of e. 
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Define also the height of the excursion e hy e = sup{e(s) : < s < C(^)}- We now obtain in the 
case y > 



(E, X E) 



rG 



{X,,Gdi),y,,ed?/,r;<Ta} 



(E, X E) 



—rG T T 
^ ''^{supft<g{?h}<a, r,&{g,g+iTaAOieg)),Xg&dv}Heg{r,-g)edy} 



(E, X E) 



dX,e 



-rGs 



{suph<s{eh} 



<a,x.edv}hv>s} n [t' < {Ta A C)(e), e{t') G dy] |^,^^_ 



Ex 



^-qs-rOsj 



(5.4) • i5;[n(r/< (T, AC)(e), e(??) G dy)] , 

where the sum in the first equality is over all left-end points of excursion intervals. The second 
equality follows from the compensation formula of excursion theory (see O Cor. IV. 11]) and the 
third from Fubini's theorem, the independence of X and i] and the fact P[r] > s] = e"''**. The 
second expectation in (|5.4|) is given by 



{5.5)E[n{r]<{TaAC){e), e{r]) € dy)] 



dy ye (0, a) 



The formula in (jS.Sp was proved in [30] (see equation (22)). 

In order to obtain the first expectation in (j5.4p first note that, since at the time LJT^ the process 
X is at its supremum level x + t (i.e. L^^ = T^^^), we must have G^-i = L~[^ P2.-a.s for all t > 0. 
Furthermore the identity X j^^i — x = t implies that if we reparametrise the integral under the first 
expectation in ()5.4p using inverse local time we find 



Ex- 



-qs-rGs J _ ^AY 



E^ 



g-$(g+r)(?j-x)]g ^ 



'I 



v — x 



dv 



V — X 



dv 



(5.6) 



The third equality follows by an application of the Esscher change of measure formula and the 
equality in (j5.6p is a consequence of the following 

{T+ < Ta} = {yt e[0,v- x] with et ^d,et< a} -a.s. 

Since Nt = #{u G [0,t] : €u > a,€u ^ d} is a Poisson process with parameter n^^'i~^^\e > a) under 
the measure Fx^'^~^^\ where n**-'''''^^ is the Ito excursion measure of X under the measure P*'-'''''^'' , 
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we obtain 
(5.7) 

(5.8) 



^_(,_,.)n*(.+^)(^>a)^ where 



$(g + r) 



VF$(g+,)(o) M^('?+0(a) 

and W$(q_|_r)(a) denotes the 0-scale function under the measure IPx ^ . The first equality in (j5.8p is 
a well-known representation of the scale function in terms of the excursion measure and the second 
equality follows from the Esscher change of measure formula W^(^q^j.^{x) = q~^^'^1'^^)w'^'^^'^\x) . 
Since Grj = t] on the set {Y^ = 0} we find in the case y = 



'rj=0,ri<Ta}' 



\ ^ {Xn>V,Yri=0,rj<Tal' 

' ^{T+<t,Yt=0,t<Ta}'^'' 



r + q 



(r + g)e-('^+'-)*/{y,=0}dt 



W 



{q+r)r 



qexjp —V' 



(a)\ _ 



where the first factor follows by combining (15. 6p . ()5.7p and (15. Sp and the second factor follows 
from |30l equation (22)]. Since l^(«+'^)(0) = W^'^\0) we thus obtain 



(5.9) 



X E)[e '^''''/|x,ed.,y,=o,r,<r4] = ^«^P 



VF(«+'-)'(a)^ 



1^('')(0). 



Ty(<?+0(a) 

Identities (^A^, ([53]), ([53]), dSS]), dSS]) and ([521) and dSj]) together imply the formula in (f3lB . 
This concludes the proof in the case of the event Aq. 

In the case of creeping, step (1) consists of the application of the strong Markov property and 
the first factor in (|3.1ip follows from an analogous calculation to (|5.ip with Aq replaced by Ac- We 
therefore assume in what follows that x > u + a. 

Step (2) does not feature in the context of Ac as there is no overshoot at t^. For the analogue of 
step (3) note that x > u + a implies (see also Figure ([TJ) the inclusion 

{Yr,=a}c{X^^>u} P,-a.s. 

The following excursion calculation, similar to (|5.4p . provides a key step: 



-qra-rG^^T 
^ ^{X^^edv,YTa=a} 



(5.10) 
(5.11) 



E, 



E^ 



{suPh<9{^4<a,^r,Gdi;}« ^' 



9>0 

CO 



-qs-rGsT _ 
^ ^{suPh<s{^h}<a,Xs£dv}^^^ 



{Ta{eg)<a^g),eg{Ta{eg))=a} 
e~''^''^'^hTaie)<ae),e{Ta{e))=a}n{de) 
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The summation in (I5.10p is over the left-end points of the excursion intervals and the equality 
in (jS.lip follows from the compensation formula of Ito's excursion theory. 

The expectation in the first factor in (jS.lip is identical to (j5.6p . To find the second factor in (jS.lip 
first note for any 2; > the following identities hold as a consequence of the compensation formula: 



E 



E 



E 



{suph<g{eh-X,,}<a, Xg<xY 



Y''Ph<L-^{'h-^h}<aj J£ 



. . dy / e ''^"+«(')/{T.+,(.)<C(.),.(T.+,(e))=a+y}n(de) 

h\<a.> Je 



The second equality follows by reparametrising the integral which arises in the compensation formula 
using inverse local time. Since Lq ^ = a.s. and the trajectories of X and are right-continuous, 
the following equality holds: 



(5.12) / e-«^"(^)/|T„(.)<c(.),.(T„(.))=a}ri(de) = lim x-iE 



al 



x„ 



This limit can also be derived from Proposition 2 in [12]. Insert identity (|2.6p into (j5.12p . apply 
L'Hopital's rule and the fact ^^(''^'(O) = ^ to obtain a formula for the second factor in (jS.lip . This 
implies (j3.1ip and therefore concludes the proof of the theorem. □ 

5.2. Proof of Theorem [S]. We start by noting that the spectral negativity of X implies the 
following identity X^^ — Xj^^ = Y^^ = a a.s. We can therefore apply the Esscher transform to (j3.20p 
to obtain 



^^qra-vG^^+uX^^j 



^-(q-^(u))Ta-rG^^j 



(5.13) E 

If a < then P \X^^ < v\ = 1 as remarked above (see also Figure [2]). Note further that 

<t-} = {r-„, <T+} for any v ^[Q,a]. 
Observe that Xrp- < v — a P-a.s. Hence the strong Markov property applied to the right-hand 

v — a 

side of (|5.13p and formula (|2.4p yield 



E" 



-ig+r-i>iu))T-_^ j- 



{TJ-_,<T+} 



E" 



~{q-^{u))Ta-rG~ 



Wi'i){a)Z^^\a) 

for all f > 0, where p = q + r — ip{u). To complete the proof we need to show that 

p^{9+0(a) 



(5.14) 



E" 



-{q-^{u))Ta~rG~ 



Wii){a)z'i'\a 



for q,r,u > 0. 
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This, together with ()3.18p - ()3.19p . imphes the formula in the theorem. 

By the compensation formula applied to the Poisson point process of excursions of Y away from 
zero, 



(5.15) E"[e^(«^'^("))^"^''&a] =E" 







e 



where n'^ is the corresponding excursion measure under the probability measure P", L a local time 
process of Y at zero and Ta and ^ are as in the proof of Theorem [T] (see page [20|) . From Bertoin [6l 
Ch. VII, Prop. 15] it follows that for any (J7)-stopping time r and ^4 S J> we have 



(5.16) 



n"(A n {r < C}) = lim VF„(x)"^E: 

x\0 



X 



An{T<T„ } 



for any u > 0. 



This limit can also be derived from [12\ Proposition 2]. We now apply (j5.16p to the stopping time 
r = T+ to obtain 



lim Wu{x)~'K 

x\0 



-{q-^{u))T+ J 



(5.17) 



^^<?-^W)(o)-i lim 



{q-il){u)) 



(x) 



for any q > ip{u). The second equality follows from (|2.3p . The third equality is a consequence of 
the well-known identity W'^^\x) = J2T=o l''W<''+^\x), which holds for all / G C, x > 0, where 
denotes the k-th. convolution of W with itself. 

The identity in (|5.15p for r = 0, the Laplace transform 
in [301 Prop. 2], and identity (|5.17p imply 

r{q-ip{u)) 



l/zi^-i'(-))^a), given 



'{a) 



for q > ip{u). This identity, together with ()5.15p and (|5.17p . yields (j5.14p for large q. Analyticity 
of expressions in (|5.14p on both sides of the identity imply (|5.14p for all parameter values. This 
concludes the proof of the theorem. □ 
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